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Abstract 

The Microwave Anisotropy Probe (MAP) is a follow-on to the 
Differential Microwave Radiometer (DMR) instrument on the Cosmic 
Background Explorer (COBE) spacecraft. The design and analysis of 
the MAP attitude control system (ACS) have been refined since work 
previously reported. The full spacecraft and instrument flexible model 
was developed in NASTRAN, and the resulting flexible modes were 
plotted and reduced with the Modal Significance Analysis Package 
(MSAP). The reduced-order model was used to perform the linear 
stability analysis for each control mode, the results of which are 
presented in this paper. Although MAP is going to a relatively 
disturbance-free Lissajous orbit around the Earth-Sun L2 Lagrange 
point, a detailed disturbance-torque analysis is required because there 
are only a small number of opportunities for momentum unloading each 
year. Environmental torques, including solar pressure at L2, and 
aerodynamic and gravity gradient during phasing-loop orbits, were 
calculated and simulated. A simple model of fuel slosh was derived to 
model its effect on the motion of the spacecraft. In addition, a thruster 
mode linear impulse controller was developed to meet the accuracy 
requirements of the phasing loop burns. A dynamic attitude error limiter 
was added to improve the performance of the ACS during large attitude 
slews. The result of this analysis is a stable ACS subsystem that meets 
all of the mission’s requirements. 
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INTRODUCTION 


The Microwave Anisotropy Probe (MAP), one of the first two Medium-Class Explorer 
(MIDEX) missions, will measure the anisotropy of the Cosmic Microwave Background (CMB), 
which is believed to be a remnant of the Big Bang, or Primordial Fireball, marking the birth of 
the universe. 1 This anisotropy was first measured by the Differential Microwave Radiometer 
(DMR) instrument on the Cosmic Background Explorer (COBE) satellite. 2 " 4 MAP has been 
designed to measure the spectrum and spatial distribution of the CMB with sensitivity 50 times 
that of the DMR and angular resolution 20 times finer, specifically 0.3° or 18 arc-minutes. These 
increases in sensitivity and resolution should enable MAP to determine the values of key 
cosmological parameters and to answer questions about the formation of structure in the early 
universe. 


MAP is scheduled to launch in the Fall of 2000 on a Delta launch vehicle, and will be 
placed in a Lissajous orbit around the Sun-Earth L 2 point using a lunar assist with phasing loops, 
reaching its final orbit approximately 100 days after launch. The MAP radiometers cover two 
fields of view (FOVs) 135° apart on the celestial sphere. To obtain a highly interconnected set of 



measurements over a large area of the celestial sphere, the MAP observatory will execute a fast 
spin (0.464 rpm) and a slower precession (0. P/sec) of its spin axis about the Sun line. The entire 
celestial sphere will be observed once every six months, or four times in the planned on-station 
mission life of two years. 

There are six ACS operational modes: Inertial, Observing, Delta V, Delta H, Sun 
Acquisition, and Safehold. Inertial mode acts as a staging mode between the other operations of 
the spacecraft; it can either hold the spacecraft in an inertially-fixed orientation or slew the 
spacecraft between two different orientations. Observing mode is used for science operations. 
Delta V mode uses the thrusters or reaction engine modules (REM) to adjust the orbit. Delta H 
mode uses the REMs to unload excess angular momentum. Sun Acquisition mode acquires and 
maintains a thermally-safe power-positive orientation of the spacecraft. Safehold mode uses 
independent attitude control electronics to put the spacecraft in a power and thermal safe attitude. 

The remainder of the paper will present various analyses of the MAP ACS. The effect of 
flexible modes, reaction wheel jitter, and fuel slosh on pointing stability performance will be 
discussed. Estimates of the environmental torques in the phasing loops, lunar swingby, and final 
L 2 orbit will be presented, along with contingency procedures for managing these torques. A 
linear impulse controller designed to improve the accuracy of MAP’s thruster firings will also be 
discussed. Lastly, a dynamic attitude error limiter was added to improve the performance of the 
ACS during large attitude slews; this will be presented as well. 

STRUCTURAL MODEL 

Reduced Order Flexible Modes 

Initially, the NASTRAN 5 model incorporated large bodies like the antenna and thermal 
reflector as point masses. The newer NASTRAN model incorporated the physical dimensions of 
these large appendages and instruments. With these new inclusions the dominant flexible mode 
frequencies increased from 1.9 Hz to 5,7 Hz. From this model a reduced order flexible model has 
been created which reduced the system to include only 3-4 dominant modes for each Spacecraft 
Mode. The associated energy for each mode was compared in MSAP 6 through the formulation of 
various summary methods: Frequency, Modal Gain, Peak Amplitude and Gregory’s Method 
which are listed below in Table 1. In addition, each mode’s singular value plots for the rigid 
body, the flexible modes, and a combination of the two were analyzed. The current estimated 
dominant mode frequencies are 5.7, 5.8, 27.7, 33.2, 36.2 and 44.9 Hz. These frequencies will 
change somewhat as the models become more precise, but the dominant modes are expected to 
remain unchanged. 

The structural system is fairly rigid and the flexible modes have only a small effect on 
the spacecraft system’s stability. As a result, the majority of analyses were done with a plant 
derived as a combination of the rigid body modes and the chosen dominant flexible modes. The 
various structural plants analyzed were: 1) the rigid body modes, 2) the rigid body modes plus 
flexible body modes, and; 3) the rigid body plus the reduced flexible body modes (lines 1, 2 and 
3, respectively, in Figure 1). The system plant was modeled as Plant = l/(Is 2 ) + Flexible modes, 
where, I is the spacecraft inertia. Since the inertia values change by less than 1% (0.1 dB) from 
the “beginning of life” to “end of life”, one stability analysis will hold for all phases of the 
mission. 
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Figure 1. Singular Value Plot of Flexible modes for MAP 




































































Stability 

Nichols stability charts were created for each spacecraft system mode in INCA 7 8 . The 
dynamic plant modeled as a structural plant multiplied by the gyro dynamics (a simple second 
order system) was included in a closed loop feedback system with a proportional derivative (PD) 
controller. By implementing the Nichols and Bode techniques with the appropriate proportional 
and rate controller gains the stability margins for all the modes were calculated. All of the modes 
except one satisfied the Guidance Navigation and Control Center design criteria of at least 12 dB 
gain margin and 30° phase margin. Typically, the gain margin fell between 14 dB and 22 dB and 
the phase margin was between 36°-81°. Since the Safehold Mode Y-axis was unable to meet the 
required bounds, it was necessary to add a low pass structural filter in that axis. The addition of a 
structural filter to the dynamic plant produced the Nichols plot shown in Figure 2. 



Figure 2. Nichols Stability Chart for Safehold Mode Y-axis 


WHEEL JITTER RESONANCE 


Since Observing mode requirements are the most stringent, it was necessary to examine 
wheel jitter resonance in that mode. This wheel jitter resonance is due to the dominant flexible 
modes resonating with the wheel imbalance torques. The forces and torques on the body due to 
the structural resonance are described as: 
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where o) r is the wheel's angular velocity, L is the length to the wheel from the gyroscope, and A sr 
and A dr are the static and dynamic imbalance of the wheels, respectively. 

MAP’s wheels have imbalance values of A sr =2.5E-5 kg-m and A dr =2.5E-5 kg-m 2 . Two 
angular velocities were considered: the first wheel frequency is near the flexible mode of 44.9 Hz 
and the second one is when the 3 wheels are running simultaneously at 25 Hz, near the 28.95 Hz 
flexible mode. The body rotation torques seen at the gyroscope due to wheel resonance were 
calculated to be 22.7-36.1 arc-sec at 44.9 Hz and between 1.04-1.65 arc-sec at 25 Hz. This 
analysis showed that the wheel jitter in resonance with the flexible mode frequencies does not 
cause problems with MAP’s Observing Mode pointing requirements. 

ENVIRONMENTAL TORQUES 

The environmental disturbance torques that may act upon the spacecraft while it is 
orbiting the Earth are as follows: aerodynamic, gravity gradient, solar radiation and magnetic. 9 
Since most of MAP’s orbit is very high, torques due to the Earth’s magnetic field should not be 
significant. 

Orbit Phasing Loops 

After MAP’s injection into orbit -perigee of 900 km and apogee of 3xl0 5 km- about the 
Earth, there are multiple phasing loops. During these loops, the thrusters will bum to increase the 
eccentricity and semimajor axis of the orbit, placing the vehicle in position for a lunar swingby. 

In this orbit, drag is the most important external torque at perigee because the density of the 
planet’s atmosphere increases exponentially as the altitude decreases. Therefore, at the perigee 
altitude, a maximum aerodynamic force and torque were calculated. In addition, the gravity 
gradient torque and resultant momentum buildup were computed for both the Earth and the 
Moon pericepses. 

Aerodynamic Torques . The drag torque and momentum build up was calculated for the portion 
of the Earth’s orbit below 1000 km altitude. Initially, the sun-side surface areas are used to 
calculate the aerodynamic forces. These areas are primarily made up of the solar panels, the 
Teflon webbing, and the middle spacecraft body hexagon. A three-dimensional model of MAP 
was created in SPAD 10 . All momentum values were calculated by multiplying the torque by the 
total elapsed time of 1 day. At the perigee altitude, the maximum aerodynamic force and torque 
were 0.0471 N and 7.68xl0' 4 N-m respectively. The orbital altitude is less than from 1000 km for 
approximately 7 minutes, which produces a momentum of 0.32 N-m-s. The sun-side SPAD 
results that were averaged during one rotation about the spin axis verified that the average torque 
is only 3.367xl0~ 4 N-m and the momentum buildup is 0.14 N-m-s. The SPAD results for the x 
and y face exposure were put into SystemBuild 11 and simulated during the portion of the orbit 
when the altitude varied from 2250 km to 300 km to 2250 km. The simulation produced data very 
similar to the previous SPAD results. The maximum x and y face aerodynamic torque was 0.012 
N-m and the momentum has a peak of 0.47 N-m-s and an average of 0.28 N-m-s. 

Gravity Gradient for Earth . The calculation of the Earth gravity gradient torque on the 
spacecraft and the momentum buildup were calculated and simulated for the portion of the Earth 
orbit below the altitude of 1000 km. The Earth’s gravity gradient momentum and an average 



torque were calculated to be 0.1744 N-m-s and 2.0x1 0" 4 N-m in the x-dircction, and 0.221 N-m-s 
and 4.1x10 s N-m in the z-direction. Gravity gradient simulation results further verified values for 
torque and momentum buildup. The Earth gravity gradient simulations displayed a maximum 
torque of 1.7xl0' 4 N-m in the x and y-axes and 1 .0x1 O' 5 N-m in the z-axis, with a maximum 
momentum of 0.15 N-m-s in the x and y axes, 0.07 N-m-s in the z-axis, and a maximum system 
momentum of 0.15 N-m-s. 

Lunar Swingby 

As previously mentioned, the phasing orbits insert the spacecraft into the proper 
trajectory for a lunar swingby. As a result, the gravity gradient torque and momentum associated 
with the Moon’s environment were calculated. These values were computed for a 20 minute 
portion of the orbit during the closest approach to the Moon, which is approximately 1000 km. 

Gravity Gradient for Moon . The torque and momentum values for the Lunar gravity gradient 
torque on the spacecraft are 3.6xlO* 5 N-m and 0.0736 N-m-s in the x-direction, and 5.0xl0* 6 N-m 
and 0.064 N-m-s in the z-direction. Gravity gradient simulation results further verified values for 
torque and momentum buildup. Lunar gravity gradient simulations displayed a maximum torque 
of 3.5xl0* 5 N-m in the x and y axes and 2.5xl0' 6 N-m in the z-axis and a maximum momentum of 
0.04 N-m-s in the x and y axes, 0.026 N-m-s in the z-axis, and a maximum system momentum of 
0.047 N-m-s. 


Libration Point 

After MAP’s lunar swingby, the spacecraft will be in orbit about the Earth-Sun L 2 point. 
The solar radiation pressure and a solar pinwheel torque associated with a solar array deployment 
misalignment have been calculated. For these computations, the primary surface areas exposed to 
the Sun are the solar array panels, the silver Teflon webbing, and a black and white painted 
hexagon. For the solar pinwheel torque, the hexagon area in the middle of the z-face was 
excluded because it will not experience an angular deflection due to deployment. 

Solar Radiation Pressure . The calculated values of solar radiation result in an instantaneous 
Solar pressure torque, T so ian and average momentum buildup, AH^ur, of 1.85xl0* 6 N-m and 0.16 
N-m-s per day. Under ideal conditions the T so iar and AH SO iar about the spin and precession axes 
average out to zero. With a misalignment tolerance of ±0.25° about the spin axes, the T so iar= 
1.07xl0' 8 N-m and AHso^ 9.3xl0 4 N-m-s for a period of 1 day. For further verification these 
calculated values were compared to SPAD generated data. The results from SPAD showed an 
average solar torque and momentum buildup per day of 1.88x1 0' 8 N-m and 0.00162 N-m-s. 

These calculations could be off by an order of magnitude because of different assumptions and 
the greater numerical precision of the SPAD method. 

Pinwheel Torque . The solar arrays may deploy improperly, resulting in a canted surface area. A 
combination of their tilted surface and the solar pressure torque could cause a pinwheel torque, 
which is a torque about the spacecraft spin (z) axis. This disturbance effect has been calculated 
for a solar array canted angle of 1° and a spin axis misalignment of ±0.5°. 

The maximum instantaneous pinwheel torque and the accumulated momentum in the z- 
axis are 1.46x1 0’ 6 N-m and 0.126 N-m-s/°-day. The average momentum buildup per day in the x 
and y-axes with a misalignment angle of 0.5° in the spin axis are 0.0736 N-m-s and 0.0645 N-m- 
s, respectively. The SPAD pinwheel torque and momentum results were produced for two cases. 



For the worst case of the sun vector parallel to the normal vector, the torque and momentum 
disturbances are 5.82x1 O' 7 N-m and 0.0503 N-m-s, respectively. For a pitch sunline angle of 
22. 5°, the torque and momentum are 5.3 1 x 1 0 7 N-m and 0.0459 N-m-s, respectively. 

The initial SPAD data were curve-fit to produce the torque and momentum equations. 
The resulting rate of momentum buildup was 0.043 N-m-s per day. Then, the values were 
modified to account for improved modeling of the MAP configuration, resulting in an increase of 
23% in the surface area facing the sun, and an increase of 1 1.7% in the moment arms. This leads 
to an angular momentum buildup of 0,059 N-m-s per day, and thus to a time of 25.5 days for the 
momentum to build up to the Observing Mode performance limit of 1.5 N-m-s. This could mean 
that a momentum dump would have to be performed every 3-4 weeks. For each simulation that 
will be run later for analysis the momentum vector in the body frame will be initialized to [0, 
±1.5, ±1.5]/V2 • The system’s momentum tolerance is 1.5 N-m-s, and the total time between 
unloading bums is required to be 90 days; the maximum solar array deployment misalignment 
angle is 0.228°. The environmental torques and momentum are summarized in Table 2. 

Table 2 


Summary of Environmental Disturbances 


Calculated Disturbances 

Toraue (N-m) 

Momentum (N-m-s) 

Aerodynamic Phasing Bum 


0.32 @ t =7min (1000 km) 

Gravity Gradient Phasing Bum 

Tx=2E-4, Tz=4.1E-5 

Hx=.174, Hz=0.221 @1 orbii 

Gravity Gradient Lunar Swingby 

Tx=3.6E-5, Tz=5E-6 

Hx=.073, Hz=0.064 @1 orbit 

Solar Pressure (instantaneous) 

Tx=1.85E-6 

Hx=0.16 @ t=1 day 

Solar Pressure (average) 

Tx=1.07E-8 

Hx=9.3E-4 @ t=1 day 

Solar Pinwheel (instantaneous) 

Tx,y,z=1 .23E-4.1 .32E-4.1 .46E-6 


Solar Pinwheel (average) 


Hx=0.073 & Hv=0.064 @1dav 

Simulated Disturbances 

Toraue (N-m) 


Aerodynamic Phasing Bum SPAD z-face 

3.147E-4 (@300km) 


Aerodynamic Phasing Bum SPAD x&y-face 

0.012 (2250km=»300krn=>2250krn) 


Gravity Gradient Phasing Bum SB 

Tx=1.7E-4, Tz=1 .OE-5 

Hx=0.15, Hz=0.07 @1 orbit 1 

Gravity Gradient Lunar Swingby SB 

Tx=3.5E-5, Tz=2.5E-6 


Solar Pressure (average) SPAD 

Tx=1 .88E-8 


Solar Pinwheel (average) SPAD 

Tx=5.31 E-7 i 


SPAD- Solar Pressure & Aerodynamic Drag 


SB- System Build 


Solution for Handling Environmental Torques 

As a response to the concern of the pinwheel torque potentially causing excessively fast 
system momentum buildup, a method was developed to unload excess system momentum while 
in Observing Mode, using a series of three “one shot” thruster firings. This algorithm is not 
currently baselined to be flown on MAP, but was developed as a contingency procedure in the 
event that system momentum build-up on station is greater than expected. The algorithm was 
designed to meet the following constraints and design goals: 

1. Unload momentum to less than 0.3 Nms while in Observing Mode. 

2. Always keep the solar arrays normal within 25° of the sunline (violations of the 
22.5±0.25° Observing Mode sunline angle constraint are permissible). 

3. Perform the entire operation during one ground pass (37 minutes) 

Using the MAP thruster pair 1 and 2 as an example, the steps in this process are as follows: 







































A. Wait until the momentum transverse to the z-axis is all in the +x axis. Fire thruster 2 to 
remove as much of this momentum as possible. 

B. After thruster firing A, wait until the sun is in the (-x,z) quadrant of the x-z plane. Fire 
thruster I or 2 (depending on the sign of the z-axis momentum) to add x momentum equal to 
the amount of momentum in the z-axis. By doing this, the system momentum vector is 
positioned such that half a precession cycle later it will be almost entirely in the -t-x-axis. 
This results in an intermediate system momentum state as much as -Jl higher than the initial 
value, but simulations show that this amount of system momentum does not pose an attitude 
control problem. 

C. After thruster firing B, wait approximately half of a precession cycle (30 minutes), until all 
of the system momentum is in the +x axis. Fire thruster 2 to remove as much of this 
momentum as possible. 

The total procedure requires a maximum of 35.5 minutes and reduces system momentum to near 
zero. The algorithm can be adjusted to use one thruster or use one of the other two thruster sets. 

FUEL SLOSH 

Propellant slosh is a mechanical effect of liquid propellants. At the moment the 
propellant mass impacts the wall of the fuel tank, it will transfer momentum to the spacecraft. If 
the tank is not at the center of mass of the vehicle, this will create an impulse moment that will 
affect the vehicle attitude. This disturbance can be significant and unacceptable, and must be 
limited through tank design. Fortunately, the MAP tank is close to the center of mass and the 
moment arm is small. In addition, a flexible membrane (diaphragm) has been placed in the tank 
to assure that the propellant remains in contact with the propellant port and to help damp out the 
sloshing effects. 

The original analysis for the ACS Critical Design Review used a model for fuel slosh 
that was found in two different reports. 12- 13 The TDRS slosh model from NASA CR 166745 
report 12 is presented in Figure 3. In that model, the slosh natural frequency, 0) n depends on the 
torsional spring constant, K, and the overall acceleration, a, of the tank: 



Figure 4 is the model used in a study for the CRAF/Cassini spherical tanks. 



Figure 3. TDRS Fuel Slosh Model Figure 4. Cassini Fuel Slosh Model 



The angular position, 6 , of the model pendulum is a representation of the position of the 
liquid bulk and the lateral offset, Zc g , of the liquid center with respect to the direction of the 
settling acceleration can be calculated as (see Figure 4): 
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Eqn. 2. Then the equations of motion for the system are defined as: 
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angular torque, T angular = JO = -J — -^-sin co n t , and 
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angular momentum, H angular = JO = — J — -^-cos 0) n t . 

From the initial analysis, the calculated fuel slosh inertia, J, is 0.7 % of the total system 
inertia and the derived natural frequency, co„ is several times larger than the system’s bandwidth 
frequency of 0.02 Hz. 


The previous reports had several questionable aspects. For example, one of the models 
did not include a stiffness factor and the equation for the natural frequency was incorrect for the 
corresponding dynamic equation. Therefore, it was necessary to derive our own scenario for each 
aspect of the fuel slosh’s movement. The five worst cases considered were: 

(1) Radial Thrust with a full tank of fuel; 

(2) Slew Maneuver with a full tank and a half full tank; 

(3) Linear Station Keeping Effects at L? (half tank); 

(4) Angular Station Keeping Effects at L 2 (half tank); and 

(5) Observing Mode Fast Spin, accelerating and constant rate (all ranges of fuel). 

Simple models were drawn for each case to obtain dynamic equations describing the motion of 
the fluid. Further research was done to find a better value for the stiffness coefficient, K and the 
damping coefficient, C. Values for K and C could be associated with MAP’s diaphragm 
thickness (0.12 in) and were extrapolated from References 14-15. The results are shown below in 
Figure 5. Calculation of the torques and momentum buildup due to the fuel slosh were done for 
each case. In addition, an angular response due to the fuel slosh torque was calculated from 
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The basic block diagram for the analytical system is shown in Figure 6. 
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Figure 6. Block Diagram to 
Analyze Angular Effect due to 
Fuel Slosh 


Figure 5. Extrapolation of MAP’s K and C with 
Diaphragm Thickness 


All the results for the various cases are presented below. 

Case 1- Radial Thrust with a full tank of fuel: 

The angular frequency was 1.05 to 1.5 Hz, the linear torque was 3.65xl0' 5 to 7.52xl0‘ 5 
N-m, the angular torque was 0.055 to 0.025 N-m, and the angular momentum was 0.004 
to 0.001 N-m-s. 

Case 2- Slew Maneuver with both a full and half full tank a fuel: 

The angular frequency was 1.04 to 1 .5 Hz, the linear torque ranged from 0.752 to 0.69 
N-m, the angular torque ranged from 4.13 to 6.36 N-m and the angular momentum 
ranged from 0.55 to 1.24 N-m-s. 

Case 3- Linear Station Keeping Effects at L 2 (half a tank): 

Radial Thrusting produces an angular frequency of 1 .64 to 2.05 Hz, a linear torque of 
8.54xl0' s to 1.13X10 -4 N-m, an angular torque of 0.011 to 0.02 N-m and an angular 
momentum ranging from 0.0007 to 0.001 N-m-s. 

Perpendicular Thrusting produces an angular frequency of 6.23 to 8.25 Hz, a linear 
torque of 1 1 N-m, an angular torque of 61 N-m and an angular momentum ranging from 
1.5 to 2 N-m-s. 

Case 4- Angular Station Keeping Effects at L 2 (half a tank): 

The angular frequency was 1.62 to 2.04 Hz, the linear torque was 0.69 N-m, the angular 
torque was between 3.76 to 4.19 N-m and the angular momentum ranged from 0.37 to 
0.52 N-m-s. 

Case 5- Observing Mode Fast Spin, accelerating and constant rate (for all ranges of fuel). 

A constant angular rate applied to the spacecraft produced an angular frequency of 0.39 
to 0.503 Hz, a linear torque between 4.6xl0' 6 -1.3xl0' 5 N-m, the angular torque was 
between 0.0014 and 0.002 N-m and the angular momentum ranged from 1.33xl0' 5 to 
1.7xl0‘ 5 N-m-s. 

An angular acceleration applied to the spacecraft produced an angular frequency of 0.39 
to 0.503 Hz, a linear torque between 0.01 and 0.017 N-m, the angular torque was between 
0.00013 and 0.00017 N-m and the angular momentum ranged from 6.4x1 0‘ 5 to 8.9x1 0‘ 5 
N-m-s. 


After this analysis, it was found that fuel slosh would not have a large affect on the 
spacecraft’s attitude. Although some of the instantaneous torque and momentum values are large, 






their overall effect was small. The natural frequency of the slosh system is several orders of 
magnitude larger than the bandwidth of the control system. The angular displacement due to 
either the linear or angular torque is from 10' 9 to 10' 5 radians. This analysis determined that fuel 
slosh should not be a concern to the MAP mission. 

CONTROLLERS 

Thruster Mode Linear Impulse Controller 

The Thruster Mode linear impulse controller was designed to improve the accuracy of 
MAP’s z-axis Delta V’s. There are some tradeoffs involved in usage of this impulse controller; 
with the impulse controller the spacecraft uses about 12% more fuel, but without it, errors during 
Delta V are about 12% larger. The Current implementation allows for a bum accuracy of less 
than 1 second, but this could be improved to 0.04 second if necessary. This excludes subsequent 
firings in the Delta H mode. 

The operational plan is to enable the impulse controller right after MAP’s lunar swingby 
during the mid-course correction (MCC) so it is available at L 2 . At L 2 , the absolute duration of 
the bums is smaller so the percentage errors tend to be higher. Without the impulse controller the 
bum is corrected with a one-sided impulse controller in both x and z. Regardless of whether or 
not the impulse controller is enabled, the current design will ensure that the bum duration will be 
at least the desired amount of time. Figures 7 and 8 shows a sample L 2 bum with and without the 
impulse controller, respectively; the commanded thruster time for each run was 60 seconds in the 
x-axis and 30 seconds in the z-axis. Notice the error caused without the impulse controller and 
the extra thruster firings to correct the error with it. 

Dynamic Attitude Error Limiter 

The dynamic attitude error limiter enables the spacecraft to meet the sunline requirement 
during a slew in Inertial Mode. When the limiter is not utilized. Inertial Mode slews that include 
high spin errors in the z-axis can cause the spacecraft to violate the 25° sun constraint, as shown 
in Figure 9. The dynamic attitude error limiter calculates an attitude error limit for each axis 
proportional to the error in that axis. This preserves the direction of the resulting Inertial Mode 
slew and prevents the spacecraft from violating the sun constraint, as seen in Figure 10. Both 
Figures 9 and 10 show a series of 45° Inertial Mode slews across the sunline with spin angles 
from 0° to 180°; with the dynamic attitude error limiter, the spacecraft does not violate the sun 
constraint for these same slews. 









Figure 9. Inertial Mode slews with no dynamic attitude error limiter 
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Figure 10. Inertial Mode slews with dynamic attitude error limiter 




CONCLUSIONS 


Analysis showed that the spacecraft is fairly rigid, and all of MAP’s operational modes 
except one satisfied the Guidance Navigation and Control Center design criteria. Since the 
Safehold Mode Y-axis was unable to meet the required margins, it was necessary to add a low 
pass structural filter to that axis. In addition, the wheel jitter in resonance with the flexible mode 
frequencies does not cause problems with MAP’s Observing Mode pointing requirements. 

Most of the environmental disturbances are manageable, but misalignment of the solar 
panels can give rise to a pinwheel torque that causes momentum buildup. This situation calls for 
a maximum allowable solar array deployment misalignment angle of 0.228°. Simulations show 
that a proposed solution of a three-shot momentum unloading in Observing Mode can be used if 
the misalignments exceed this limit. 

The natural frequency of the fuel slosh is several orders of magnitude larger than the 
bandwidth of the controller. Therefore, the angular displacement due to the fuel slosh’s linear 
and angular torque is small, and fuel slosh should not be a concern to the MAP mission. 

The thruster impulse controller and the dynamics attitude limiter were added to improve 
the performance of the attitude control system. The end result of the analysis is a MAP controller 
that meets attitude control system requirements with margin. 
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